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The role of the algebraic method has long been understood in shedding light on the topological 
structure of sets. However, when the set is a simplicial complex and host to a dynamical process, 
in particular the trajectory of a canonically distributed system in thermal equilibrium with a heat 
bath, the algebra re-enters. Via a theorem of Levy and Dynkin, there is a correspondence between 
a system's energy process at equilibrium and a sequence of n— chains on the state space. 



FUNDAMENTALS AND REVIEW 
Thermostatics 

Let the set {u>i, u>2, . . . , wjv} be the state space 
for a random process {xt ■ t € [0, oo)}. Accord- 
ing to the results of yj, under the assumption 
that the "position" process Xt is a trajectory of 
a canonically distributed subsystem in equilib- 
rium with a heat bath at temperature 9, there 
is a well defined energy process 

h{ Xt )=6\0g(^-r) (1) 

V p(x t ) ) 
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where, Atk is the characteristic time associated 
with the k th state during a cycle and 

j_ 

LT = (Aii At 2 • •• At N ^ " . 

For a simple example see 3] . Sample Isotherms 
as seen from the time domain are shown in figure 
□ 

With these microscopic quantities in hand, 
one is in a position to calculate the free energy 
and many other important quantities. The in- 
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n=[p lP 2---p N ) (2) 

and for k = 1,2, ...,N 

Pk = p(x t 6 OJk) ■ 

Further Q, if the N— state process Xt has a 
characteristic Carlson depth (in the finite state 
Markoff case this is simply the L 1 norm of the 
cycle between visits to a rare state) we have 

(3) 

I 

ternal energy 

U = p • h 

will be of particular importance in the next sub- 
section. 

Thermodynamics 

It may happen that there are changes to the 
observed probabilities, changes in the length of 
the characteristic time scale or both, if it is as- 
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FIG. 1: Surfaces of Constant Temperature for a Two 
State System 



sumed that these changes are occurring in a sub- 
system that is in a state of continual equilibrium 
with a heat bath (or perhaps the observer shifts 
her attention from one subsystem to another all 
in established equilibrium with the same bath), 
the notions of work done, heat gained or lost by 
the system become applicable. 

In equilibrium these changes are given by the 
product rule for differentiation (the First Law) 
as 

U= p^h + p^h (4) 

work rate heat rate 

In the nonstationary case equation @ trans- 
forms in the obvious way from II into II(i). Sim- 
ilarly h(xt) becomes h(xt,t), etc. The calcula- 
tion of the rates involve time derivatives of these 
quantities and are easily computed. The vectors 
appearing in the dot products on the RHS of 
equation J3J are typically high dimensional. 

Fortunately, in the statistical mechanical con- 
text, the work rate admits a low dimensional 
description as well. Typically this is system de- 
pendent. Famous examples include pressure and 
volume changes or magnetization and changes in 
an applied field. The low-dimensional descrip- 
tion of the heat rate is often system indepen- 
dent and is given in terms of temperature and 
entropy. 



Frequently, as the microscopic dynamics often 
occur on state spaces with a topology these same 
macroscopic phenomena may also admit a sec- 
ond, low-dimensional interpretation in terms of 
their interplay with features of the state space. 
The purpose of this paper is to illustrate the 
twofold nature of the thermodynamics. The 
connection is accomplished via the subsystem's 
trajectory through its path space and the path 
decomposition of Levy and Dynkin. 

Dynkin-Levy Decomposition of the Path 

This subsection focuses some elementary as- 
pects of the martingale theory on applications 
to the energy process. In particular, the context 
is that of a canonically distributed subsystem's 
trajectory through its state space. For readers 
unfamiliar with the martingale point of view an 
excellent reference for this section is Q|. 

Recall that according to the Levy Formula, 
see for example @, the martingale decomposi- 
ton of a Markoff process is given by: 

g(x t ) - g(x ) - \ ^2 9x 3 k (g(k) - g{x s )) ds 

J ° k^ xs 

here 

P(x s + =k\x s ) 
q X3 k= hm 

s+^s S + — S 

is the generator of the process. 

Of particular interest in the sequel is the case 
when g(-) is the energy function h(xt, t) and the 
temperature is held constant. Under these con- 
ditions the Dynkin-Levy kernel becomes 

. (. n(«+) , n(«) \ 

Note: it may also be that the generator is itself 
a function of time. 

On the Choice of State Space 

As the main purpose of this paper is to bring 
out the connection between the physical and and 
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algebraic descriptions of the subsystem's trajec- 
tory it is worthwhile to choose a state space that 
is meaningful in either context. 

Exchange rules in the context of urn models, 
, 0, U, are frequently used as state spaces 
to study diffusion, spin flip dynamics, granular 
materials, phase transitions, etc. Further the re- 
sulting state spaces are often simplexes or sim- 
plicial complexes which have long been used as 
building blocks in the study of algebraic geom- 
etry and topology. 

A particularly simple member of this family of 
exchange models is the Gordon-Newell queue . 
For an interpretation of that process in terms of 
its thermodynamics see Q. Figure |21 illustrates 
the algebraic topology associated with its dy- 
namics. Despite the potential for an arbitrarily 
large number of non-degenerate (energy levels 
are distinct) microstates, both the thermody- 
namic and algebraic descriptions of the energy 
process are three dimensional. The connection 
will be more fully developed in the sequel. 

• 

• • • 
• • • • 



FIG. 2: A basic example of a 2— simplex derived 
from the state space for Gordon-Newell exchange. 
The 0— chain, 1— chain and 2— chain shown (color 
coded) are derived from the martingale decomposi- 
tion of the local energy process. 



Elements of the Gordon-Newell Process 

The state space of a Gordon-Newell process 
may be represented as a matrix. The dimen- 
sion of the column space is the number of urns. 



The rows represent the states accessible to the 
system and may be used as state labels. 

For an instantiation based on M balls and m 
urns, the cardinality of the state space is simply 
the number of ways to put M balls (indistin- 
guishable) in m urns. In equation (JSJ an ex- 
ample of the state space matrix and row label- 
ing scheme are presented. The left hand side of 
equation (J5J defines the list of state labels as- 
sociated with a space based on three urns and 
three balls (a total of ten states). The right hand 
side of the equation lists the states themselves. 
The dimension of the column space is equal to 
the number of urns. 

In the sequel there will be essentially two 
types of vector dot products, those describing 
macrostructure and those describing microstruc- 
ture. For clarity, dot products with respect to 
the low dimensional row space will be denoted 
with a circled circle notation, for example (see 
equation J5J) 

a(7) © r\ = lrji + 2t] 2 + 0r) 3 . 

Dot products in the column space and other high 
dimensional vectors will use the standard " ■ " 
notation. 

The labels a(l),a(4) and a(10) are the 
0— simplexes. These are the states with all the 
balls in one urn or, alternatively, the states with 
two urns empty. The states a(2) and £*(3) be- 
long to the same 1— simplex, namely the set of 
states with the first urn empty and the others 
non-empty. The label a(6) state belongs to the 
interior of the region. That is, the set of states 
with no empty urns. 
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A basic version of the exchange rule consists of 
an to x to global routing table (a transition prob- 
ability matrix for an embedded discrete process) 
and to— timescales that characterize (globally) 
the rate of activity of the individual urns. 

Let Pij denote the probability that the dis- 
crete process transitions a single ball from urn 
"i" to urn "j". The stationary vector (eigen- 
value 1) of the transition matrix will be denoted 
{711, 7T2, . . . , 7r m }. The urn rates will be denoted 
{qi,q 2 , ■ ■ -,q m }- 

Any component of a stationary measure of the 
process may be given in terms of these parame- 
ters: 

( — ) Ql ( — ) Q2 ■ ■ • ( %m \ a "> 



qi 
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where, {oc\, 012, ■ ■ • , a m } are entries in that row 
of the state space matrix associated with the 
state label a(-). 

A slightly more flexible version of the ex- 
change rule allows the routing table and rates 
to vary locally from state to state. In addition, 
the parameters may be functions of time. 



THE CONNECTION WITH THE HEAT 
RATE PROCESS 



of the exchange rule parameters and was intro- 
duced in the last section. The middle expression 
follows using the substitution 



77. = -log(— ) 

q. 



The bottom expression gives a stationary mea- 
sure in terms of energy and temperature. 

As is apparent from equation J3J , the temper- 
ature may vary (via dilatation of the time scales) 
even though the state probabilities themselves 
are constant in time. When both the tempera- 
ture and probabilities are constant in time 



h(x t ,t) 



ft (af t ) 

Qi] ® (a(x t ) - a) 



(7) 



where a is the arithmetic mean of the rows of 
the state space matrix. For the example given 
in equation JSJ this vector is {1,1,1}. 

Dynkin-Levy kernel for the energy process is 
given by 



The Constant Temperature, Time 
Stationary Equilibrium 

With the state energies and temperature 
given in terms of the canonically distributed 
subsystem's statistics and timescales by equa- 
tions CJ and a stationary measure /z(-) of 
the exchange process may be expressed (compo- 
nentwise) in a few ways 



^ qi ' ^ qi ' ^ <?m ' 



fj,(a) = < 



g— (7)1 ai +7)20:2 +—+Vm°lm) 



(6) 
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Recall that the top expression on the RHS of 
equation © is the stationary measure in terms 



1} ® (a(x s ) - a(k)) 

k^x s 



h(a(x B ))-h(a(k)) 



In what follows the Dynkin-Levy coefficient at a 
site "*" will be denoted dl(*). 



Physical Significance of the Dynkin-Levy 
Kernel for the Energy Process 



Under the hypothesis of this subsection, con- 
stant temperature, time stationary equilibrium, 
the Dynkin-Levy kernel takes the form given in 
equation |JSJ. It may be shown that this corre- 
sponds to the heat process associated with the 
energy process h(xt). For 
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(dl(x s ))ds = y^y^ p(x s ,s) P Xsk (s,s + ) 

x s k 

x [a(k) - a(x s )] ©r) 

= J] p(k, a+) [a(k) — a] © q 
fc 

= (p(m, s+) — s)) 
x [a(m) — a] ® jj . 

So that, 

(dl(x s )) = Y^p(m, s)h(m, s) . (9) 

m 

This is of course the heat rate's contribution to 
the "First Law", equation |0J. 

Algebraic Significance of the Dynkin-Levy 
Kernel for the Energy Process 

Inspection of equation (JSJ shows that the 
Dynkin-Levy coefficient at a site depends on the 
a differences with the states to which the pro- 
cess may transition, weighted by the probabil- 
ity rate that the transition takes place during a 
small time s — > s + . 

Recall from the previous section that in the 
most basic version of the exchange rule there is 
one global matrix (site independent) that gov- 
erns the transitions of the embedded discrete 
process and one set of inverse characteristic 
times that governs the rates of the urns. 

A typical local neighborhood of an interior 
point "A" with state label a(A) = {a, b, c} is 
shown in figure The quantities relevant to 
the calculation of the Dynkin-Levy coefficient 
at site "A" are listed in the table. 

Notice that the table entries are invariant with 
respect to the actual location of site "A" . That 
is, the Dynkin-Levy coeficients are constant on 
the interior of the simplex. A similar situa- 
tion holds in the case "A" belongs to one of the 
boundaries. 



TABLE I: A break down of the summands that com- 
prise a Dynkin-Levy coefficient: possible transition 
states, a— index differences and corresponding small 
time probabilities of transition. See figure [3] 
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a(k) - a(A) 
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1,0,-1 


q-iP-ii ds 
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0,1,-1 


93^32 ds 


D 


-1,1,0 


qiPi2 ds 


E 


-1,0,1 


giPi3 ds 


F 


0,-1,1 


q 2 P23 ds 


G 


1,-1,0 


q 2 P2i ds 



(a+1,b-1,c) (a+1,b,c-1) 




(a-1,b,c+1) (a-1,b+1,c) 

FIG. 3: Typical Neighborhood of an Interior Point 
in the State Space 



Consider the 1— simplex of states that have 
the first urn empty and the other two urns non- 
empty (for a picture see the orange 1-boundary 
in figure |2J). The typical neighborhood and ta- 
ble for this boundary are obtained from those 
shown above for the interior points by simply 
deleting points "D" and "E". The same table 
and figure result for all points on this face. In 
this way a real valued coefficient is associated 
with each 1— simplex in the complex. The set of 
these coefficients specifies a particular 1— chain. 

The reader should be able to convince her- 
self that under the conditions of the basic ex- 
change rule: one global (site independent) tran- 
sition matrix and one global set of characteristic 
times for the urns, the Dynkin-Levy coefficients 
for the energy process specify n~ chains in higher 
dimensions as well. 
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THE CONNECTION WITH THE WORK 
RATE PROCESS 

It may happen, see equation© or figure 2] 
that the probabilities vary in time but the sys- 
tem remains on a constant temperature surface. 
In this situation the quantities n_ , q_ and r\_ from 
equation © become functions of time. The 
Dynkin-Levy kernel of the energy process takes 
the form (see the appendix, equation l(T%|) ) 

dl(x s ) = - 2J Qx s k(s) 

k^x s 

x (h(a(x s ), s) - h(a(k),s + )) 

= 9 T)(s) ® (a(x s ) - a) 
+9 ^2 1x s k{s) 

k^x s 

xij(s) ® (ot(k) - a(x s )) 

(10) 

Comparison of ifTUI) with equation (JSJ reveals 
an additional term 

9fi(s) © (a(x s ) - a) . 

This site dependent term is clearly the local 
work rate at x s . For (see equation 

6J2 P(x.) [r{(s) © ( a (x s ) -&)] 
= 9 ^p(m) [r}(s) © (a(m) — a)] 

m 

= P h. (11) 

SUMMARY 

Under the condition of global (site indepen- 
dent) exchange rules on a simplicial complex, 
a correspondence has been shown between the 
chain groups (with real coefficients) of the state 
space and the thermodynamics. In the case 
of constant temperature, time stationary equi- 
librium, the Dynkin-Levy kernel ties the sys- 
tem's local heat rate to specific members of those 
groups. 



Under the same conditions, in the case of con- 
stant temperature but time varying dynamics, it 
is shown that, due to the presence of a site de- 
pendent (local) work rate term, the correspon- 
dence with the chain groups is lost. Thus, the 
connection has been made between transitions 
from one thermal equilibrium to another and 
from one sequence of p-chains on the state space 
to another. 



APPENDIX A: DERIVATION OF THE 
DYNKIN-LEVY KERNEL FOR THE 
ENERGY PROCESS 

Let h(xt,t) denote the energy of the state xt 
at the time t. That is, the Hamiltonian for the 
system may have an explicit time dependence. 
The collection of all system trajectories (the po- 
sition processes {x t , t g [0, oo)}) will be denoted 
f2. Let V denote the probability on path space 
the describes the trajectories of the position pro- 
cess. 

The energy difference between the an "initial" 
state Xt at time to and a "final" state xt m at 
time tmi may be decomposed into a telescoping 
sum: 

Shll™ = H x t m ,t m ) ~ h(x to ,t ) 

= [h(x tm ,t m ) - h(x tm _ 1 ,t m - 1 )] 

+ [H x t m -i,t m -i) - h(x tm _ 2 ,t m -2)] 
+ ... 

+ [h{x tl ,t 1 )-h(x t0 ,t )} . (12) 

As is customary, let T t be the a— algebra gen- 
erated by paths in f2 up to time t. For simplic- 
ity, is supposed that To contains the informa- 
tion about the explicit time dependence of the 
Hamiltonian. 

Also, as is customary, let E{*\Tt) denote 
the projection operator (conditional expecta- 
tion) that projects its argument onto the set of 
square integrable functions on the measurable 
space (Q, Tt) using the measure V. The accu- 
mulated difference between the actual trajectory 
and its projections may be approximated as fol- 
lows 
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-L h \to = h ( x t m ,t m ) - h(xt m _ i ,t m _ 1 ) 
- E(h(x tm ,t m ) 

+ h ( x t m _ 1 ,t m _ 1 ) - h(xt m _ 2 ,t m _ 2 ) 
-h{xt ,t m -*)\Ft_S) 



+ [h(x tl ,tt) - h(x to ,t )] 

- E(h{x tl ,tx) - h(xt ,t )\T to ) (13) 

It is a standard result that _L h \\ m is a martin- 
gale with respect to the filtration J- t previously 
described. 

In the case that the position process is 
Markoff, the small time contributions to the ker- 
nel may be computed as 



= [h(x tm , t m +i) - h(x tm , t m )} 
x(l - 0(t m+1 - t m )) 

+ ^2 [h(k,t m +x) - h{x tm ,t m )] 

X E tm fe(tm-1, t m ) (14) 

more precisely: 

E ( h(x tm+1 ,t m+1 ) - /i(x tm ,i m )|^ r tm ) 
_ H x t m , t m+ i) - h(x tm , t m ) 

, , \t>m-\-l v<m ) 

+ ^ [ft.(fc,t m+ i) - ^(x tm ,t m )] 

X Q.a: tTO fc( (Wl-*m) ■ (15) 

Assuming that 

h(k,a + ) - h(k,s) 

is small, for example h as a differentiable func- 
tion of time, this is 



E ( h(xt m+1 ,t m +i) ~ K x t m ,tin)\Ft m ) 
_ h(x tm ,t m+ i) - h(x tm ,t m ) 



2J [h(k,t rn ) - h(xt m ,t r 



(16) 



Under these conditions the accumulated af- 
fects over the the time interval [0, t] are 



' ( d h(x s ,s) 
~d~s 



At constant 0, in terms of the exponents from 
equation Q 

9 I r)(s) © — a] 



(18) 



x q (s) ds 



At constant 9, in terms of fancier expo- 
nents resulting from the application of time 
varying microforces (local exchange rules and 
timescales) 



ds I r](s) © [ct{x s ) — a] 



+ S(x a , s) © a(x s ) - S(*, s) © a(*) 



17(a) © [ot(k) - a(x s )] 



+S{k, s) © a(k) — 5(x s , s) © a(x s ) 



(19) 
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